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Phonon self-energy corrections have mostly been studied theoretically and experimentally for 
phonon modes with zone-center (q = 0) wave-vectors. Here, gate-modulated Raman scattering 
is used to study phonons of a single layer of graphene (1LG) in the frequency range from 2350 
to 2750 cm -1 , which shows the G* and the G'-band features originating from a double-resonant 
Raman process with ij / 0. The observed phonon renormalization effects are different from what 
is observed for the zone-center q — case. To explain our experimental findings, we explored 
the phonon self-energy for the phonons with non-zero wave-vectors (q ^ 0) in 1LG in which the 
frequencies and decay widths are expected to behave oppositely to the behavior observed in the 
corresponding zone-center q — processes. Within this framework, we resolve the identification 
of the phonon modes contributing to the G* Raman feature at 2450 cm -1 to include the iTO+LA 
combination modes with q ^ and the 2iTO overtone modes with q — 0, showing both to be 
associated with wave-vectors near the high symmetry point K in the Brillouin zone. 

PACS numbers: 73.20.Hb, 73.22.-f, 78.30.Na, 78.67.Ch 
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Elcctron-phonon (el-ph) interactions are responsible 
for many important effects in condensed matter physics 
[l|, |2fl. In particular, the phonon self-energy, which is 
mainly due to the el-ph coupling, is a remarkable effect 
which contributes to both the phonon frequency and de- 
cay width renormalizations due to the creation (annihi- 
lation) of electron-hole (e-h) pairs through phonon ab- 
sorption (emission). For non-zero Fermi energy values 
(Ep ^ 0), phonon self-energy corrections are needed to 
explain a set of well-known effects, such as the Kohn 



anomaly [U 



the Peierls transition [l|, ls|, LZ|] , polaron 
formation [1|, |8|, |9j , and other types of phonon renormal- 
izations and perturbations caused by e-h creation (anni- 
hilation) due to phonon absorption (emission) 10l4l5|. 




Single-layer graphene (1LG) has linear electronic en- 
ergy dispersions E(k) around the non-equivalent high 
symmetry points K and K' in the Brillouin zone as a so- 
lution of the Dirac equation which gives massless particle 
behavior around K (K 7 ) [g, |lg, [l7 1 . However, one cannot 
properly solve the electronic and vibrational structure for 
most carbon materials near the Dirac points when consid- 
ering the adiabatic approximation, which disregards the 
ionic motion of the carbon ions lg, Il8l421|. When the 
adiabatic approximation cannot be applied [iHj], the el- 
ph interactions must take into account non-adiabatic pro- 
cesses, which give rise to important and strong phonon 
self-energy corrections [lj|. Within second-order pertur- 
bation theory, the phonon self-energy can be approxi- 



mately described as [1| 



nK,Ef) = 2^ 



IVkk'l 2 



kk' 



hbJa — E, 



kk' 



*7q/2 



X(fh-fe) (1) 



where k and k' are, respectively, wave- vectors for the 
initial and final electronic states; q = k — k' is the 
phonon wave-vector; E^, = (££, — E^.) is the e-h pair 
energy; cj q is the phonon frequency; 7 q is the phonon 
decay width; fh(fe) is the Fermi distribution function 
for holes (electrons) and Vkk' gives the el-ph coupling 
matrix element information. For a specific u; q , the 
phonon energy correction due to its self-energy is given 
by fiw q — fto; q = Tie[Tl(u!q, Ep)], which is the real part 
of EqJTJ where Huj® is the phonon energy in the adia- 
batic approximation. Likewise, the decay width 7 q is 
given by the imaginary part Im[n(w q , Ep)] of Eq.[TJ Note 
that while 7 q /2, which is self-consistently determined 
by Im[II(u; q , Ef)] |lll-ll4l|. is always a positive quantity, 



the quantity hw q 



K 



Rc[Il(w q , Ep)] is positive if 



hujq > E^_, or negative if ftw q < E^.,. If foj q - fiw|! > 0, 
a phonon hardening in w q occurs while if fiw q — ftw q < 0, 
a softening occurs. 

These phonon renormalizations occur any time we have 
non-zero matrix elements Vkk' and occupied (unoccu- 
pied) initial (final) states, in the sense that an electron- 
hole pair can be created (annihilated) by a phonon ab- 
sorption (emission) process as a perturbation. As defined 
in Fig.[T] there are two types of electron-phonon interac- 
tions, namely intra-valley (AV) (Figs.QJa) and (cX) and 
inter- valley (EV) (Figs.QJb) and (d)) processes [22|. For 
an AV process, the initial and final states both occur 



within the region close to a K [K'] point, while for inter- 
valley processes, K is connected to K' (or K' to K), 
respectively, by a q ^ phonon. Thus the AV (EV) pro- 
cess corresponds to T (K) point phonons. The phonon 
wavevector q for an AV (EV) process is measured from 
the r (K) points and can assume both the q — or q / 
conditions. Because of momentum conservation require- 
ments, the q = condition connects electronic states 
with the same wave- vector fc, while the q ^ condition 
connects electronic states with different k values. 

In the past, most of the discussions of the phonon 
self-energy renormalizations have been for zone-center 
phonons (r-point) with q = 0, which can be appreciated 
by observing the G-band Raman feature evolution in 1LG 
as the Fermi level is varied (see Fig.[3|)0, [l(| [HI, UM, [23 - 
251 ] . In the present work, we use gate-modulated reso- 



nant Raman spectroscopy (RRS) to address the effect of 
a gate voltage (Vg) on the phonon self-energy (EqfT|), for 
1LG systems, in cases where ? / (AV and EV pro- 
cesses). These are cases which have not been sufficiently 
studied previously. Here, we study the double resonance 
Raman frequency ranges between 2350 and 2850 cm -1 , 
which contains the G* and the G'-band spectral features 
as shown in Fig.^a) 0, 0, Ej]]. We show that the 
phonon renormalization for the q + 1 and q = phonons 
measured from the K (K') points in the Brillouin zone 
gives a different Fermi-energy dependence to that for the 
q = r-point phonons. We apply these differences in 
behavior to show that the G* feature is composed of two 
Raman peaks which behave differently from one another 
as \Vg\ is varied. 

The graphene flakes used in our experiments were ob- 
tained by micro-mechanical exfoliation of graphite over 
silicon substrates with a 300 nm thick layer of Si02. 
Next, e-beam lithography was performed to pattern our 
devices. Then, thermal evaporation of Cr/Au (5 and 
80 nm, respectively) was done. For each Vg value, 
RRS spectra were taken with a 532 nm wavelength laser 
source in the backscattering geometry using a 100X ob- 
jective. The laser power measured from the objective 
was 1.5 mW. The spectra were analyzed using Lorentizian 
line-shapes from which frequencies and decay widths were 
extracted. Figures[2ja) and (b) and Figs.[3][a)-(e) show 
the experimental results. Note that RRS provides infor- 
mation about both the electronic and vibrational struc- 
tures, while the Vg variation allows for control of Ep. 
Due to the difference in behavior between the q = and 
q =/= processes, this combination of techniques provides 
a precise way to verify the assignments of either over- 
tones and/or a combination of phonon modes. It is worth 
commenting that the results presented here are observed 
for several samples and for several locations on the same 
sample. 

The G* and G' features were intentionally chosen for 
this discussion because: (1) they are the most prominent 
double-resonance Raman features in the graphene spec- 



trum, offering a convenient platform, together with the 
G-band feature, to observe experimentally the two differ- 
ent types of phonon renormalizations, one found for the 
q = phonons and the other for q ^ 0, and (2) by using 
these different phonon renormalization effects, we have 
solved a long-time discussion in the literature about the 
G* feature, showing that the G* feature is composed of 
both the iTO+LA (q = 2k) and 2iTO (q = 0) modes, 
both measured from the K-point. In the literature, the 
G* feature around 2450 cm -1 has been assigned to ei- 
ther the iTO+LA phonon combination mode (q = 2k 
EV process) [271 ]. or the 2iTO phonon overtone mode 
(q = EV process) |26j, awaiting a more definitive as- 
signment. Note the possibility of two types of double 
resonance conditions, q = and q = 2k, for a phonon 
with momentum q and an electron with momentum k 
[H 0]. We show here that both the 2iTO (q = 0) 
and the iTO+LA (q = 2k) modes are Raman active and 
therefore both are expected to appear. The G' (or 2D) 
feature at 2670 cm -1 is widely known to be an overtone 
of the iTO phonon mode (q = 2k) 0, S 113] ■ It gives a 
dispersive phonon frequency as a function of laser energy 
^Lascr which exhibits the value of 103cm -1 /eV [33j. The 
iTO+LA (q = 2k) combination mode presents a disper- 
sion of -(16±l)cm -1 /eV (measured in this work), while 
the 2iTO (q — 0) overtone mode (also measured in this 
work) does not disperse (see Fig.^b)). FigureEJa) shows 
that indeed the G* feature is asymmetric, suggesting that 
it consists of two Lorentizians peaks rather than just one. 

EquationQ] has previously been explored for the cases 
where the phonon momentum q vanishes (q — 0) for the 
AV intra-valley process. In these cases, ftu; q — fiu' = 
Re[n(w q , Ep)} and 7 q = Im [II ( w q , Ep)] give the results 
illustrated for the renormalization of u; q and 7 q , respec- 
tively, in the insets of Figs.[2ja) and (b), respectively. 
Observing these insets, we conclude that the renormal- 
ization is a minimum when Ep — 0. When \Ep\ becomes 
approximately half of the phonon energy (E q /2 = hujq/2) 
and greater, the renormalization becomes strong because 
of decreases in the energy denominator in Eq.[T] [14} . As 
an example of phonon renormalization when q = for 
the AV process (see Fig.Qla)), the ujg and jg variations 
of the G-band Raman feature are shown in Figs. [21(a) 
and (b), respectively, as Vg is varied. The G-band fea- 
ture corresponds to the first order q = iTO (in-plane 
transverse optical) and iLO (in-plane longitudinal opti- 
cal) phonon branches around the T point. As regards 
the log an d 7g behaviors as a function of Vg , the experi- 
mental results (Figs. [Ha) and (bV) are in good agreement 
with literature predictions [4j, [5j, I23I425J , which show a log 
hardening and jq narrowing when Vq increases. Next, 
we report the new experimental results for phonons cor- 
responding to the cases q = EV (inter- valley) and q / 
AV/EV. 

Both the G* iTO+LA mode at 2450-53 cm -1 and the 
G' mode at 2670-73 cm -1 are EV double-resonance Ra- 
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FIG. 1. (Color online) (a) possible (Ef = 0) and not possible (Ef 7^ 0) AV processes for e-h pair creation (annihilation) due to 
phonon (with energy hwq) absorption (emission), (b) possible (Ef = and -Ef 7^ 0) q — (measured from the K point) EV 
processes, (c) not possible (Ef = 0) and possible (Ef 7^ 0) AV processes for electron-hole pair creation (annihilation) due to 
phonon (with energy ftu q ) absorption (emission) when the phonon wave- vector is not zero (q =/= 0). (d) possible EV processes for 
electron-hole pair creation (annihilation) due to phonon (with energy hojq) absorption (emission) when the phonon wave- vector 
is not zero (q ^ 0). 
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FIG. 2. (Color online) (a) and (b) show, respectively, the fre- 
quency lug hardening and decay width 70 narrowing for the 
G-band Raman feature as a function of gate voltage. The 
insets in (a) and (b) are theoretical pictures of the Ef depen- 
dence of Re[n(u; q , Ef)] and Im[n(u; q , Ef)] for an AV q = 
process. The Ef values on the upper scales are obtained with 
Ef=!i\vf I ^TrCgVc/e, where C g =ll5aF/(im 2 is the gate ca- 
pacitance, e is the electron charge and |«f| = 1.1x10 6 m/s is 
the Fermi velocity. The Re[II(a; q , -Ef)] and Im[n(a; q , Ef)] 
values in the insets in (a) and (b) were normalized to their 
maximum values to illustrate the concept of o; q hardening and 
7 q narrowing. 



man processes with q 7^ (see Fig.[3Ja)) and, as shown 
in Figs.[3][c)-(e), they show a different behavior when Vq 
increases compared to the behavior observed for the AV 
q = process. Starting with the G'-band feature (2iTO 
EV process with q = 2k measured from the K-point), 
it is seen that its frequency ujqi decreases with increas- 
ing \Vg\ (Fig.|3fc)), while its decay width 7c increases 
with increasing \Vg\ (see the inset in Fig.[3Jc)). Here, we 
see that the same behavior is observed for the iTO+LA 
mode frequency cjjto+la (EV process with q = 2k mea- 
sured from the K-point), as shown in Fig.[3][d), and for its 
decay width 7iTO+LA, as shown in Fig.(3][e). The 2iTO 
G* feature at 2470-73 cm" 1 , which is a q — EV pro- 
cess around the K point (see Fig.QJb)), is observed in 
Figs.[3]Jd) and (e) to show a frequency W21TO & n d decay 



width 72iTO that almost do not change with increasing 
Vg\- This behavior shows that the 2iTO (q = 0) mode 
14J couples weakly to the electronic states in graphene 
and therefore its phonon self-energy corrections are neg- 
ligible. These findings show a different w q and 7 q depen- 
dence with \Ep\ when compared to the AV q = phonon 
process (see Fig. [2]). 

To explain our experimental findings, a phenomeno- 
logical formulation for the phonon self-energy for the 
EV <z = and AV/EV (j ^ processes in single-layer 
graphene is presented as follows. In the case of AV pro- 
cesses for the q = phonons (Fig.QJa)), which applies for 
the G-band feature (see Fig.[3fa) and (a)), the renormal- 
izations are a minimum when E-p — 0, which implies 
a phonon frequency (w q ) softening and phonon decay 
width (7) broadening, and both w q hardening and 7 q 
narrowing with increasing \E-p\, as shown in the insets in 
Figs.[DJa) and (b). More precisely, the phonon renormal- 
ization shows its smallest values when Ep = ±hujq/2, 
which represent two singularities in Eq.[TJ as shown by 
the black solid curve in the inset of Fig.^a). These 
singularities give rise to what is commonly known as 
Kohn- anomalies. However, it is difficult to experimen- 
tally observe these Kohn-anomalies near room temper- 
ature because of thermal excitations (relaxations). The 
absence of such divergencies is also attributed to a large 
non-uniformity in the density of carriers due to foreign 
molecules and charge traps in the substrate 25J. There- 
fore, the condition \Ep\ » hojq/2 usually applies and 
the gray solid line-shape in the inset of Fig.^a) is ex- 
pected to be observed, which means that fiw q — hcu® is 
always positive when \Ep\ increases [13J, 1 141] . Indeed, 
for q = the E^., must be larger than 2|£ , p|, and if 
2\E-p\ >> hojq, no real e-h pair due to phonon absorp- 
tion(emission) will be created and hio^ — ftw q is positive. 
For simplicity, the discussion hereafter is based on the 
assumption that \Ep\ » ftw q /2 for non-zero \Ep\. 
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FIG. 3. (Color online) (a) The experimental G* and the G' bands as they appear in the resonant Raman spectrum. The 
asymmetric G* feature is a combination of the iTO+LA (q = 2k read from the K point) mode and the 2iTO (q — read from 
the K point) mode. The G' mode is an overtone of the iTO mode (g = 2k). For illustrative purposes, the signal of the G* 
feature was multiplied by a factor of 10 and the Lorentizian profiles used to fit the spectrum are shown in constructing (a), 
(b) The frequency dispersion of the G* peaks as a function of laser energy (i?Lascr) and shows that the iTO+LA (q — 2k) is 
a dispersive mode [-(16 ± l)cm _1 /eV], while the 2iTO (q = 0) is non-dispersive [-(0.4 ± 0.3)cm _1 /eV] [2^, Hg. (c) The 
gate voltage Va dependence of the 2iTO (q — 2k) oj g i and "/ G i (inset in (c)). (d) and (e) show, respectively, the w q and 
7q dependencies on \Ef\ seen for the iTO+LA and 2iTO modes, (f) shows illustrative predictions for the Vb-dependence of 



the phonon frequency correction 



(black solid line) and the corresponding decay width 7 q (grey dashed line) when 



q t^ 0, both as a function of Ef. In (c), (d) and (e), the Ef values were obtained as shown in Fig. [5] The Re[II(Wq, Ef)] and 
Im[n(o; q , Ef)] values in (f ) were normalized to the maximum value to illustrate the concept of aj q softening and 7 q broadening. 



Note that, a different behavior is expected for the 
q = phonon (measured from K-point) in the EV process 
shown in Fig.QJb), which explains the G* 2iTO mode be- 
havior as \Vg\ is varied. The probability that an electron- 
hole pair exists at Ep — (upper line of Fig.fTJ) is small 
since the density of states at Ep almost vanishes and 
this implies that no £j?-dependence for the frequency w q 
and width 7 q is expected. When |-Ep| increases (lower 
line in Fig.Q]), the probability for a q — phonon to 
connect equal k and k' states (k = k') increases, which 
means that the number of e-h pairs increases and thus the 
renormalization effect could be appreciable with increas- 
ing \Ev\- However, for q — EV processes, the E^, will 
be always approximately null (E^., ~ 0). This means 



that 



will be a small correction and, therefore, 



small w q and 7 q renormalizations are expected for any 
|.Ef| value (weak ^-dependence). 

By considering phonon modes with q ^ (AV and 
EV processes) as shown in Figs.QJc) an d (d), the phonon 
wave- vectors are either around the T point or around the 
K point. This case explains the G' and the G* iTO+LA 
mode behaviors as \Vg\ is varied. Since the phonon en- 
ergy dispersion for graphene has a much smaller slope 
than that for the electronic energy dispersions 17[ , there u; q — w q and 7 q must behave oppositely to the behavior 
is essentially no coupling between 9 ^ phonons with observed in the insets of Fig.[5J This opposite behavior is 



e-h pairs if Ep = and therefore phonon renormaliza- 
tions in this case do not take place in a resonant way, 
i.e. where E^i — E^ = hujq. This means that no phonons 
with 9/O can connect electronic states with different 
k and k' at Ep — and as a result, the matrix ele- 
ments Vkk' in Eq.Q] are close to zero and essentially no 
self-energy corrections occur. Precisely speaking, in the 
case of the EV process (Fig.QJd)), the scattering is pos- 
sible for Ep — 0, but the density of states (DOS) is very 
small. However, when Ep ^ 0, phonon modes with q =/= 
can now connect electronic states with different k and 
k', giving rise to a strong electron-phonon coupling and 
then the phonon self-energy corrections become impor- 
tant (see Figs.QJc) an d (d))- Such a q ^ phonon will 
be observed in the double resonance Raman spectra [28j , 
and in this case one can observe defect-induced features, 
overtones or combination modes [a, lla, LL2| ■ 

For g / AV and EV processes, the requirement 

> 2\Ep\ no longer exists and for most of the sit- 

> huiq will be satisfied and 

will be negative. Now for the q ^ AV 

and EV cases, only phonon softening and decay width 

broadening are predicted to occur. This requires that 



rpeh 

uations the condition E^ 
then 



illustrated in Fig.|3ff), where it is seen that the frequency 04147. 
correction w q — w" (black solid line) must become more 
negative with increasing \Ep\ while the decay width (grey 
dashed line) must increase with increasing \Ep\. 

In summary, the widely studied intra- valley AV q = 
case [4J, [5j, I23I425J shows that when Ep = the phonon 
renormalization is a minimum. The phonon frequency 
then shows a hardening while the decay width shows a 
narrowing when \Ep\ increases. Here, we have studied 
the phonon self-energy correction for phonon modes with 
q / (AV and inter-valley EV processes) and q — (EV 
process) from a theoretical and experimental point of 
view. In the q ^ cases (oppositely to what is observed 
for the q — AV process), the phonon renormalization 
is a maximum when E-p = and the phonon frequency 
softens while the decay width broadens with increasing 
|_Ep|. Namely, while the decay width 7 q is always posi- 



tive, the frequency correction w q 



is positive for the 



<7 = AV case but negative for the q ^ AV/EV cases. 
For the EV q = case, Efy ~ and a weak and small w q 
and 7 q dependence with Ep is expected. In this context, 
gate-modulated resonant Raman spectroscopy of over- 
tones and of a combination of phonon modes provides a 
powerful technique to assign the phonons participating 
in the formation of each overtone or combination mode, 
to identify whether a Raman feature is associated with 
the q — or q ^ processes and to determine how a 
given phonon mode is coupled to the electronic states 
of single-layer graphene. As shown in Figs.|3](a)-(f), we 
applied these combined techniques to study the G* and 
G' modes, which are the most prominent double reso- 
nance q ^ Raman features in the graphene spectrum. 
Our theoretical approach satisfactorily explains the ex- 
perimental results and within this framework, we also 
showed that the G* mode is an asymmetric peak com- 
posed by both , the iTO+LA combination mode, which 
is an EV q = 2k process with aw q -w?<0 renormaliza- 
tion and the 2iTO overtone mode, which is an EV q = 
process with a weak phonon renormalization, thereby re- 
solving a long-time discussion in the literature. 
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